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ABSTRACT 



The X-ray flux of Nova V2491 Cyg reached a maximum some forty days after optical maximum. The X-ray spectrum at that time, 
obtained with the RGS of XMM-Newton, shows deep, blue-shifted absorption by ions of a wide range of ionization. We show that the 
deep absorption lines of the X-ray spectrum at maximum, and nine days later, are well described by the following phenomenological 
model with emission from a central blackbody and from a collisionally ionized plasma (CIE). The blackbody spectrum (BB) is 
absorbed by three main highly-ionized expanding shells; the CIE and BB are absorbed by cold circumstellar and interstellar matter 
that includes dust. The outflow density does not decrease monotonically with distance. The abundances of the shells indicate that they 
were ejected from an O-Ne white dwarf. We show that the variations on time scales of hours in the X-ray spectrum are caused by 
a combination of variation in the central source and in the column density of the ionized shells. Our phenomenological model gives 
the best description so far of the supersoft X-ray spectrum of nova V2491 Cyg, but underpredicts, by a large factor, the optical and 
ultraviolet flux. The X-ray part of the spectrum must originate from a very different layer in the expanding envelope, presumably 
much closer to the white dwarf than the layers responsible for the optical/ultraviolet spectrum. This is confirmed by absence of any 
correlation between the X-ray and UV/optical observed fluxes. 

Key words, stars: individual: V2491 Cyg - Novae - X-ray binaries - cataclysmic variables - Interstellar Medium - chemical abun- 
dances - interstellar dust - X-ray spectroscopy 



1. Introduction 

The sensitivity to soft X-rays of the EXOSAT satellite enabled 
the first disco very of X-rays from a classical nova, GQMus, 
in April 1984 (Oe gelman etailll984l) . The ROSAT satellite de- 
tected it almost eight years later, in Feb 1992, and found it to 
have a very soft black body-like spectrum with kT e ff ~ 28 eV 
(Oeg elman et aljfl993h . ROSAT observations of other novae in 
our galaxy show ed that only a small fraction of novae are de- 
tected in X-rays (lOegelman et alJIl 993b lOrio et al.ll200ll) . 

The spectrum of GQ Mus is very soft, comparable to the 
spectra of a new class of supersoft X-ray sources (SSS) discov- 
ered with ROSAT in the Large Magellanic Cloud, with tempera- 
tures < 100 eV in bl ackbody fits to their spectra (Triimpe r et al.l 
I199U iGreineretal 1 ll991l) . The favored explanation for the X- 
ray spectra of the supersoft sources, as for GQ Mus, is hydrogen 
burning at the surface of a white dwarf, of material accreted from 
a companion star. The class of supersoft X-ray sources shows 
a large variety, some sources being permanent and others tran- 
sient, and with different types of c ompanion stars (see review by 
Kahabka & van den Heuvel lT997). Observation campaigns with 
ROSAT, XMM-Newton and Chandra in M 31 show that most 
of the supersoft so urces are classical novae (iPietsch et aLl l2005: 
iHenze et alj2010h . Recentlv. lHenze et alj d2011l) argued that that 
the high fraction of novae without detected SSS emission might 
be explained by the incomplete coverage. 

The model fits of the supersoft X-ray spectra have in- 
creased in sophistication over the years. Black body fits in- 



dicated highly super-Edd ington luminosities for the ROSAT 
sources. iHeise et all (1 19941) showed that the fit of local thermo- 
dynamic equilibrium (LTE) model atmospheres to the ROSAT 
data brought their bolometric luminosity below the Eddington 
limit, for temperatures < 60 eV. For temperatures > lOOeV, 
Hart mann & Heisel (fl997h showed that the differences between 
bolometric luminosities for ROSAT supersoft sources derived 
from blackbody, LTE and NLTE models are much smaller, the 
main difference being that NLTE models may give apprecia- 
ble flux above 0.4 keV (i.e. in the 'h ard' ROSAT band) . A new 
level of sophistication is reached by iPetz et all (120051) who fit 
NLTE model spectra of expanding atmospheres, thus taking into 
account the outflow inherent in the nova phenomenon. Even 
the most sophisticated model fails to fit the observed contin- 
uum of the nova V4743 Sgr at A < 30 A, and als o does not 
match the depth of the abso rption lines (Fig. 3 of IPetz et ail 
120051: iRauch & Werner! I2010I) . With regard to the X-ray spec- 
tra of V2491 Cyg, the best published atmosphere model does 
not reproduce the observed line velocities and does not repro- 
duce the observed line depths, as clearly visible in Fig. 11 of 
Ness et al. (201 1), for example in the N vn ls-2p line at 24.8 A. 
Atmosphere models so far also fail to reproduce the wide range 
of ionization levels simultaneously present in the observed spec- 
tra: e.g. O v to O vm in V2491 Cyg (see Figs.GBbelow). 

What is the reason for this failure? If there are shocks in the 
expanding nova envelope, the assumption of a standard atmo- 
sphere that density decreases monotonically outwards, breaks 
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down. The rapid X-ray variability observed in some novae in- 
dicates that the expanding envelope and/or the luminosity from 
the deeper layers are not stationary and/or not spherically sym- 
metric, invalidating another assumption in atmosphere models. 
In this paper we therefore investigate a pilot model of a rather 
different nature, in which separate expanding shells absorb the 
spectrum of an underlying central source. In this first study, we 
describe the central source as a blackbody spectrum. We apply 
our simple model to the XMM-Newton RGS spectra obtained 
from the nova V2491Cyg to investigate whether this model 
gives a better description of the continuum and line absorption; 
whether the short-term variability is due to changes in the source 
or in the absorption; and whether the model allows determination 
of the chemical abundances in the expanding envelope. 

Nova V2491 Cyg was d iscovered on April 10, 2008, at 
V = 7.7 dNakano et al.ll2008l) . The optical flux declined from 
this peak in the course of the following months, interrupted by 
a brief second maximum at the end of April. Ultraviolet and X- 
ray fluxes of the V2491 Cyg were obtained with the Swift satel- 
lite: the ultraviolet flux declined in tandem with the optical flux 
(measured through May and June), whereas the X-ray luminos- 
ity rose from a marginal detection on April 1 1 to a peak more 
than a thousand times brighter some 40 days later (Page et al. 
1201 Oal) . Li ghtcurves of optica l, ultraviolet and X-ray are shown 
in Fig. 4 of iPage et all (1201 Oal) . In the Swift X-ray range (0.3-8 
keV) the spectrum peaks strongly near 0.4-0.5 keV during the 
X-ray maximum and has a flatter (i.e. harder) energ y distribu- 
tion b efore and after X-ray maximum (see Fig. 3 of Page et aD 
[2010ah . 

From a relation between maximum magnitude of novae and 
their rate of de cline, a distance of 1 0.5 kpc has been estimated for 
V2491 Cyg bv lHelton et all (|2008). The optical flux was already 
declining during the first measurements, and therefore the actual 
flux maximum may have been higher than the flux of the first 
measurement: this would imply a smaller distance. 

The interstellar reddening E(B - V) = 0.43 has been es- 
timat ed from the ratio of the O I lines at 0.84 and 1.13 mi- 
cron dRudv et al.ll2008l) . and may be converted to a hydrogen 
column Nh — 2.4 x 10 2 5 m 2 using the relation determined by 
iPredehl & Schmittl (Il995l) . 

The X-ray spectra described in this paper were obtained with 
XMM-Newton just before the X-ray maximum, and some nine 
days later, when the X-ray flux had declined by an order of mag- 
nitude. 

Our paper is structured as follows. In Section 2 we describe 
the XMM-Newton data and provide a brief description of the 
four spectra that we will analyze. In Section 3 we describe a ba- 
sic model and apply it to the spectra. In Section 4 we discuss 
several variant models and their applications. In Section 5 we 
discuss our model fits and in Section 6 we summarize our con- 
clusions and outline prospects for improvements to the models. 

2. Data and brief description of the spectra 

The data discussed in this paper were obtaine d during two obser - 
vations with XMM-Newton, as described in INess et al.l (1201 lh . 
The first observation started at day 39.86, measured from the 
initial optical discovery on 2008 April 10.73 UT, and lasted 
39280 s. During this first observation the RGS countrate varied 
strongly: from about 15 cts/s during the first 7500 s it dropped 
to a minimum near 3 cts/s lasting some 6000 s, and then rose 
to a maximum of ab out 18 cts/s for the final 16000 s of the ob- 
servation. We follow INess et aTl (1201 ll) in defining three spec- 
tra for this first observation, Spectrum A. 1 obtained during the 



initial plateau; Spectrum A. 2 during the minimum flux; and 
Spectrum A.3 d uring the final high-flux plateau, (see Fig. 3 of 
INess et al]|201 lh . The second observation started on day 49.62, 
and showed a countrate which varied only mildly around 3 cts/s 
during the 31700 s observation. All data from this observation 
are collected in the fourth spectrum, Spectrum B. 

We analyze the 7-38 A (1.77-0. 33keV) first or der spectra of 
the R GS detector, using SPEX version 2.02.03 dKaastra et al.l 
fl996l) . 

We show the four spectra derived from the XM M-Newton 
RGS observations in Figs. Q]|4] (see also Figs. 6,7 of INess et al.l 
2011). The most prominent features are the neutral oxygen ab- 
sorption edge near 23 A with the corresponding ls-2p absorp- 
tion line at 23.5 A. The neutral nitrogen ls-2p absorption line 
near 3 1.2 A is also seen, but the neutral iron L2 and L3 edges 
that one expects near 17.15 and 17.45 A are weak. Ness et al. 
(2011) identify ls-3p and ls-4p lines of Ovu near these wave- 
lengths. Indeed, absorption lines from O v to O vm are detected, 
blueshifted with respect to the rest frame wavelength; examples 
include O vm near 1 9 A, O vn at 2 1 .5 A, and N vn at 24.5 A. The 
blueshift indicates that the absorbing matter is expanding, and 
the wide range of ionization levels indicates that several expan- 
sion shells are involved. 

Lines of highly ionized neon and oxygen are also detected at 
the rest frame wavelengths; examples are Ne x at 12.1 A, Ne ix 
at 13.4A, Ovmat 19. OA, andOvnat21.6 A. A strong feature is 
also seen near 29.3 A, close to the forbidden N vi line at 29.5 A. 



3. The basic model 

As a first attempt to describe the spectrum F a we use a multi- 
component model which may be written symbolically as 



Fi = 



\ 



Fa,bb Y\ Xi + Fa,c 
i=1 



X g -X d (1) 
The components in this model and their parameters are 



- the blackbody flux F a,bb(R, ^eff) characterized by the radius 
R and effective temperature T e g, 

- three photo-ionized absorbers X,(A^,,o, V;, cr,-, (A / (9),), each 
characterized by the column density for oxygen A^q, ion- 
ization parameter outflow velocity Vj, velocity width cr,, 
and the abundances of various elements with respect to oxy- 
gen (A/0)j. The number of absorbers is in principle free: we 
add absorbers until the fit does not significantly improve any- 
more (as indicated by the Ax 2 ', see TableQ}. In this case, this 
leads to at least three absorbers. 

- the flux from a collisionally ionized plasma 
F^rXEM, T,o- c ,(A Ne )c,(MglNe) c ), characterized by 
the emission measure EM = J n e rixdV, the temperature T, 
the turbulent velocity crc, the neon abundance (A^)c, and 
the abundance ratio of magnesium to neon {Mg/Ne)c, 

- the absorption by interstellar or circumstellar gas 
X g (Nf{,T g ,A g ), characterized by the hydrogen column 
Nh, the temperature T g , and the abundances for each 
element A g , and 

- the absorption by interstellar dust Xd{Nd,o), characterized by 
the column of the oxygen in the dust Nd.o- 

To correctly describe the interstellar absorption X g with high 
spectral resolution, we use the SPEX model hot, putting its tem- 
perature T g at a low value (near 1 eV), suitable for the (nearly) 
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Fig. 5. The unabsorbed SED used for V2491 Cyg (solid line and 
data points). We also show the unabsorbed blackbody emission 
(blue dashed-dotted line) and the observed absorbed fluxes (red 
dotted line). 



neutral component of the ISM. This model describes the trans- 
mission through a layer of collisionally ionized plasma. We give 
abundances of the ISM normalized on proto-solar values; here 
and in the foll owing, when we refer to p roto-solar values we use 
those given byEodders & Palmel (l2009l) . 

The detailed form of the Fe L edge near 17.2-17.7 A and the 
OK edge near 22. 7-23. 2A may be aff ected by the presen ce of 
these elements in molecules in dust (Pae rels et al.ll200 H). The 
absorption Xd by s uch molecules is described in SPEX with amol 
dPintoet alj|2010l) . In the spectra of V2491 Cyg, the iron feature 
is weak, and thus we only constrain oxygen compounds in the 
interstellar dust. 

A highly ionized absorber X, is modeled in SPEX with the 
xabs model, which gives the transmission of a slab of photo- 
ionized material. Before this model is applied, one must com- 
pute the ionization levels, and for this we use vers ion C08.00 of 
Cloud y (for a description of an earlier version, see iFerland et alj 
1 19981) . Cloudy computes the ionization balance of a gas for the 
ionization parameter £ = L/(nr 2 ), a measure of the number of 
photons per particle (for a light source with luminosity L at dis- 
tance r of the gas with number density n), and a given spec- 
tral energy distribution (SED). The form of the energy distribu- 
tion is important, because a soft irradiating spectrum leads to a 
different ionization equilibrium than a hard irradiating spectrum 
with the same luminosity. The energy distribution that we use for 
V2491 Cyg is shown in Fig. [5] It has been computed through all 
the a vailable arc h ival d ata from IR to X-ray energies: IR fluxe s 
from iNaik et all (l2009h optical from lHachisu & Katol ((2009), 
UV and so f t X-ray from lNess et al.l (1201 ll) . and hard X-rav from 
Page et aD (l2010bl) . We have taken their fluxes and models in 
each wavelength band and corrected for the Galactic interstellar 
extinction by assuming a Nu = 2.5 x 10 25 irT 2 , which represents 
an average between our fits and those found in the literature. 
The total luminosity of this SED is 1.15 x 10 33 W, of which 
the X-rays, roughly described as a black body with the tem- 
perature of HOeV (blue dash-dotted line in Fig. [5), contributes 
6.23xl0 32 W. 

The abundances of the highly ionized absorber obviously 
also are important for the calculation of the ionization structure. 
In the computation of the ionization levels we use proto-solar 



abundances. In principle, once we have determined the abun- 
dances of the highly ionized absorbers from our spectral fit, the 
computation of the ionization structure should be iterated for 
these new abundances. In this pilot study, however, we have not 
performed this iteration. 

To limit the number of free parameters, we limit the num- 
ber of independent abundances in the basic model as follows. 
The flux F,\,c of collisionally ionized plasma dominates only at 
wavelengths A < 15 A in the spectra of V2491 Cyg and at these 
wavelengths only magnesium and neon contribute to the emis- 
sion. We determine the Mg and Ne abundances from this part of 
the spectrum, and assume that the abundances thus determined 
for the collisionally ionized plasma are also valid for the highly 
ionized absorbers Xj. Conversely, we determine the column of 
O, and the number ratios of C, N, Si, S, Ar, Ca, Fe with respect 
to oxygen from the absorption lines of the highly ionized ab- 
sorbers X/, assuming that these ratios are the same for all three 
highly ionized absorbing shells. We then assume that the ratios 
found for the highly ionized absorbing shells also apply to the 
collisionally ionized plasma. 

The reason for listing abundances with respect to oxygen is 
that these are much better constrained by the observed spectrum 
than abundances with respect to hydrogen, because hydrogen 
only contributes to the X-ray spectrum as a source of free elec- 
trons. The ionized hydrogen of the hot nova shell only affects the 
spectral normalization but not the shape in the RGS spectra. 

3.1. Application to Spectrum A.3 

We apply the basic model first to the spectrum with the highest 
signal-to-noise, SpectrumA.3. In order to allow the use of x 1 
statistics, we bin the RGS data in bins of about 1/3 of the spectral 
resolution of ~0.065 A of the first order RGS spectrum. To ensure 
a minimum of 10 counts/bin we in addition rebin the raw data at 
the shortest wavelengths. This leads to rebinning by factors 10, 
5 and 2 in the ranges 7-1 1 A, 1 1-15 A, and 15-38 A, respectively. 

In a first step we fit a preliminary model, in which we ignore 
the highly ionized absorption shells X, and the dust absorption 
lines Xj, This provides us with initial estimates for the black- 
body and collisionally ionized plasma, which we then enter as 
starting values in the full model fitting. The value that we derive 
for the temperature of the collisionally ionized plasma is deter- 
mined mostly from the spectrum at A < 15 A, and more specif- 
ically from the ratio between the Neix and Nex ls-2p lines in 
this region. For this reason we fix the plasma temperature T in 
the full model fitting at the value found from the preliminary fit. 

The resulting values of the parameters for the full model are 
listed in TableQ] and the resulting model spectrum is shown in 
Fig.E] 

The radius of the blackbody emitter, 17600 km, is sev- 
eral times the radius of a white dwarf (which decreases from 
^9000 km for a 0.6 Mp white dwa rf to ^3000 km at 1.3 M H ; see 
Eq.(27) of lNauenbergl(ll972l) . and lProvencal et al.l (120021) 1. This 
implies that the photospheric surface lies in the expanding lay- 
ers. With the temperature of 1.4 x 10 6 K, the blackbody has a 
bolometric luminosity Lbb = 8.5 x 10 32 W. 

To describe the high-ionization absorption lines, we require 
three distinct X, components. An outflow with constant veloc- 
ity and constant mass flux M, has nr 2 cc M/v and therefore £ is 
constant with radius. If the velocity increases with radius, as in a 
homologous expansion, we expect ^ to increase with radius. This 
agrees with our finding that the shell with the highest ionization 
parameter has the highest velocity. However, our finding that 
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the shell with the intermediate velocity has the lowest is puz- 
zling. We think that this result may be due to systematic errors 
in the velocities of weak shells 2 and 3, caused by the difficulty 
in separating the velocity dispersion and shift in the absorption 
lines of these relatively weak shells (see also Sect. 13.1.21 . 

The lines of C vi, Ar xvi and S xiv are weak relative to the 
oxygen lines, and this results in number ratios C/O, Ar/O and 
S/O much lower than solar. For Ca we can only give an upper 
limit to the Ca/O ratio. We estimate the hydrogen column by 
assuming a proto-solar O/H ratio. 

The relative strengths of the Ovny plus FeiL edge at 17.2- 
17.7 A and the oxygen features at 22.7-23.2 A indicate an Fe/O 
ratio lower than proto-solar. This suggests that Fe may be de- 
pleted by dust. Indeed, we prefer to describe the broad absorp- 
tion feature at 22.7-23.3 A with a dusty-molecule component of 
silicates plus a gaseous component of On, with about 10% of 
the total oxygen in each of these two components. The alterna- 
tive explanation with low-ionization O ii-iv in gaseous form only 
is less appealing as it requires a very high velocity dispersion. 



3.1.1. Deviations 

There are some features that limit the quality of the fit. Most 
of them show broad emission-like features and are much more 
relevant in the low continuum spectra. Indeed, in the high flux 
spectra these features are well hidden by both strong continuum 
and absorption lines. The strongest lines are at 18.4A, 19.0 A, 
19.2 A, 21.5 A, 21.9 A, 26.5 A, 27.5 A and 29.35 A. Most of the 
emission lines are blue-shifted less than the absorption line of 
the same transition, thus the emission should originate from a 
broader region of the shell. This would explain the whole line 
shape similar to a P Cygni profile. There are also a few absorp- 
tion lines not reproduced by the model at 1 9 . 8 A, 25 . 8 A, 29 . 1 A 
and 2 9.7 A, which have already been reported by Nes s et aD 
(1201 lh as unidentified lines. These lines could represent an ad- 
ditional slower (-2 000 km s~') expanding layer as their pro- 
jected wavelengths are consistent with O v, Ca xm, Ar xm and 
Fe xvi/Alxi transitions. However, as mentioned above, these ab- 
sorption lines are so much weaker than the strong ones, that the 
inclusion of an additional layer does not provide a significant 
improvement to the fit. 



3.1 .2. Systematic errors 

In Table [T] only the relatively small statistical parameter errors 
are given which are likely much smaller than the combined sta- 
tistical plus systematic uncertainties. The statistical errors are a 
proxy for the quality data, while systematic uncertainties are the 
main limitation for the accuracy of our results. A few param- 
eters are highly model-dependent or have a degenerate value. 
Both the velocity dispersion and the blue-shift of the shell lay- 
ers have this problem. The complexity of the absorption lines 
plus the presence of weak emission lines complicate the model- 
ing as the different components could mix with each other. This 
can happen to the components that refer to the weak unresolved 
emission lines (layer 2 and 3). Thus, the conclusion of different 
expansion velocities in two low-ionization layers is less certain 
than, e.g., the fact that the highest-ionization layer has a higher 
velocity shift and broadening. 



3.2. Application to Spectra A.1, A.2 and B 

In fitting the other X-ray spectra of V2491Cyg, which have 
rather lower signal to noise, we fix the velocity dispersion of 
the collisionally ionized plasma, as well as the (relative) abun- 
dances of this plasma, of the highly ionized absorbers X,-, and of 
the cold absorption component X g to the values determined for 
Spectrum A. 3. The resulting values of the other parameters are 
collected in TableQ] and the spectral fits are shown in Figs.[T|[2] 
and[4] respectively. 

The blackbody radius of Spectrum A. 1 is similar to that of 
Spectrum A. 3, but the temperature is lower by a quarter, and as 
a result the bolometric luminosity is lower by 50%; within the 
RGS band the blackbody flux is lower by 60%. For the collision- 
ally ionized plasma, the temperature is the same in the spectrum 
Spectrum A. 1 as in Spectrum A. 3, but the emission measure, and 
therefore also the luminosity, is lower by 45% in the former. 
As regards the highly ionized absorption shells, all parameters 
are similar in the Spectra A. 1 and A. 3, with the exception of the 
column of the shell with the highest ionization, 1.5 times big- 
ger in A.l, and the column of the shell with the middle ioniza- 
tion, which in Spectrum A.l was 0.5 times that of A. 3. Finally, 
the columns of the cold absorbers X g and Xj both are markedly 
higher, by factors 1.3 and 1.2 respectively, in Spectrum A.l with 
respect to A. 3. 

The blackbody radius of Spectrum A.2 is only 10000 km, 
and thus closer to the white dwarf surface. One might expect 
a higher temperature there, but the observation shows a lower 
temperature, ten percent lower even than the low temperature 
of Spectrum A.l. The luminosity of the blackbody accordingly 
is just 10% in SpectrumA.2 of that in SpectrumA.l. The emis- 
sion measure of the collisionally ionized plasma in Spectrum A.2 
is intermediate between those before in SpectrumA.l and after 
in SpectrumA.3, at unchanged temperature. The ionization pa- 
rameter of all three high-ionization absorption shells is lower 
in Spectrum A.2. The column of the shell with the highest 
ionization has increased by a factor of almost three between 
SpectrumA.l and SpectrumA.2; in contrast the columns of the 
other two shells have plunged. The column of the cold absorb- 
ing gas is the same in Spectrum A.2 and SpectrumA.l, but the 
column of the dust much higher in Spectrum A.2. 

The large changes in the spectral components well within the 
eleven hours of the first XMM-Newton observation are striking; 
we will investigate these some more in Sects. |4]and|5] 

Almost ten days later, during the second observation of 
V2491Cyg with XMM-Newton, the radius of the blackbody 
component has shrunk to 1900 km. This remarkably low value 
would imply that only part of the surface of the white dwarf 
has this high temperature; or alternatively that the white dwarf 
is more massive than 1.3 M . The luminosity of the blackbody 
component has declined to Lbb = 3.8 x 10 30 W. The emission 
measure of the collisionally ionized plasma has decreased, but 
its temperature has increased, and as a result its luminosity now 
is Lcie = 1.45 x 10 28 W, to a large extent within the RGS band 
(see Table[T). 

The ionization levels of the three highly ionized absorbers 
have dropped, as have the columns of shells 2 and 3, between 
Spectra A.l or A. 3 and Spectrum B; the column of the highest 
ionization shell, shell 1, has increased. The column of the cold 
gas absorber has dropped by some 10%, but the column of the 
dust component has increased by a factor 1.7. In general, the 
relative changes of the parameters of Spectrum B with respect 
to those of Spectra A.l and A. 3 are similar qualitatively, if not 
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Table 1. RGS best fit parameters to the four observations (see also Fig . [T] El [3] and©. 
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<a) The CIE temperature is kept frozen to the value estimated with a local fit to the 7 — 15 A spectral range. {b> We can only provide upper-limit for 
Ca relative abundance. w All the errors shown are statistical, systematic effects are not considered here (for further information see Sect l3. 1.2t . 

4. 1 . The hot absorption shells and the bolometric luminosity 

In the basic model the hot ionized absorber consists of three 
shells, all with the same abundances, but with distinct ionization 
parameters 

As an alternative we have fitted a continuous absorption 
measure distribution (AMD), using the SPEX model warm, to 
Spectrum A. 3 (see Tableland Fig. [6). We find that an initially 
continuous distribution in the course of the fitting procedure col- 
lapses to discrete components, compatible with the findings of 
our basic model. The model confirms the presence of at least two 
discrete ionization ranges with discrete velocities, represented 
by log £ = 3.50 - 3.75 and log £ > 4.75, which are consistent 
with what we have found with separate xabs components (see 



quantitatively, as the changes of the spectrum with the lowest 
countrate during the first observation, i.e. Spectrum A.2. 



4. Variant models 

To test the robustness of some of the results from the basic mod- 
els, we have fitted some variant models. The topics that we wish 
to investigate in particular are the nature of the hot absorption 
shells and the cause of the changes during the first observation. 
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Table 2. Absorption measure distribution. 
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' The values refer to the warm model fits to spectrum A.3. 




Fig. 6. Absorption measure distribution (see also Table |2). 



Table[TJi. The warm model uses the same velocity for the whole 

range of but the observed spectrum shows a relation between 

ionization and velocity, and this causes the warm model fit to be 

decidedly worse than the basic model Qf 2 /dof = 17300/ 1480). 
We have also fitted Spectrum A. 3 with some models in which 

we fit the abundances separately for each shell; this does not 

improve the fits significantly. 

While the limited number of absorbers suggests the presence 
of discrete shells with different velocities and a non-monotonous 
£ distribution, we note that any number of absorbers is possible 
and therefore a continuous distribution. The two strong peaks in 
the absorption measure distribution suggests the presence of two 
main absorbers. However, the high velocities make difficult to 
believe that the shells might be stable throughout the first month 
after the outburst and argue in favor of a continuous distribution. 

The temperatures of the hot absorbing shells imply that hy- 
drogen is fully ionized. This means that the hydrogen column of 
these shells is badly constrained, and for this reason we give the 
abundances of the various elements with respect to oxygen, fix- 
ing the hydrogen column to the proto-solar hydrogen-to-oxygen 
ratio in the basic model. The main effect of the presence of hy- 
drogen is that photons that initially move from the central source 
in our direction are scattered into another direction. In the fully 
spherically symmetric case, these photons are replaced by pho- 
tons moving initially in other directions that are scattered into 
our direction by Thompson scattering. The slab model that we 
use for the absorbing shell assumes that each photon from the 
central source that is scattered within the shell does not reach 
the observer, and is not replaced by photons from initially other 



for the X-ray spectrum of Nova V2491 Cygni 



Table 3. Radii and luminosities of the central blackbody emitter. 







flfkm) 




Lbb (10 30 


W) 


r eff (ev) 




basic 


noH 


noH 
no He 


basic 


noH 


noH 
no He 




A.1 


21000 


7000 


5600 


440 


47 


31 


91 


A.2 


10000 


4800 


4100 


57 


13 


9 


81 


A.3 


17600 


3780 


2790 


849 


39.2 


21.4 


121 


B 


1900 


1500 


1300 


3.8 


2.4 


2.2 


95 



The values refer to the fits to the four spectra, for the basic model, and 
for variant models in which the hydrogen content, or the hydrogen and 
helium content of the highly ionized absorption shells are set to zero. 
Note that the parameters for the basic model are correct only for a highly 
non-symmetric case; in the spherically symmetric case the correct val- 
ues for the basic model are those listed above for the no H, no He case. 



directions. In the presence of hydrogen, the number of electrons 
in the shell is much higher than in the absence of hydrogen, and 
a much higher number of photons is scattered. For a given ob- 
served number of photons, the slab model with hydrogen leads 
to a much higher number of source photons incident on the shell 
than the model without hydrogen. 

Whether this implementation of the model is valid in the case 
of V2491 Cyg is not clear. The rapid variations between spectra 
A.l, A.2 and A.3 may suggest that spherical symmetry does not 
apply. If spherical symmetry does apply, however, the number 
of photons should not be corrected for scattering, i.e. the lumi- 
nosity of the central source will be much lower than found with 
the slab model. To investigate the magnitude of this effect we 
consider first a variant model in which we set the hydrogen con- 
tent of the three shells to zero, but keep the helium. To good 
approximation, the electron scattering of the photons from the 
central source does not depend on energy within the RGS range, 
and thus the parameters of the fit are not affected, other than that 
the luminosity of the central source is reduced. In particular, the 
temperature of the blackbody is the same as in the basic model. 
Thus we compute the luminosity for the blackbody emitter for 
the model in which the photon number is not corrected for elec- 
tron scattering in the absorption shells by fixing all parameters to 
the values of the basic model, and then recompute the luminosity 
of the emitter for spectrum A.3. This luminosity, at 3.92 10 31 W, 
and with it the radius, at 3780 km, are much smaller than in the 
basic model. 

This reasoning applies equally to a second variant model 
where we set the hydrogen and helium content both to zero. For 
spectrum A.3 the temperature again is the same as in the basic 
model, but the luminosity and radius drop to 2.1410 3I W and 
2790 km, respectively. Table[3] lists the variation of the parame- 
ters with assumed hydrogen and/or helium content for the slab 
model, or alternatively the variation with increasing spherical 
symmetry. 

In computing the ionization equilibria with Cloudy we as- 
sume that the ultraviolet/optical flux is emitted by the central 
source. If instead it is emitted elsewhere, it may not contribute 
to the ionization of the shells. For this reason, we have calcu- 
lated the ionization balance of the shells when the UV fluxes in 
the adopted SED were decreased by up to two orders of mag- 
nitude. The bolometric luminosity decreased by less than 5%, 
while the newly-fitted ionization parameters of the shells are sys- 
tematically lower by 15-25%. The non-monotonous ^ structure 
as seen in Fig. [6] is still preserved and the other main physical 
parameters like column densities and abundances are consistent 
within the errors. The overall scenario thus is unaffected. 
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4.2. The changes between Spectra A.2 and A.3 



Table 4. Derived mass loss rates and shell parameters. 



The dramatic changes in brightness of the X-ray source during 
the first observation with XMM-Newton are explained in the ba- 
sic model by a drop in the luminosity of the central blackbody 
emitter, combined with changes in the column of the ionized ab- 
sorption shells. The luminosity drop of the central source is the 
result of a smaller radius and a lower temperature (see Table[3]l. 
In the variant models discussed in the previous section, where 
the hydrogen and helium contents of the absorption shells are set 
to zero, the drop in luminosity between A.3 and A.2 are much 
less pronounced, since in these models the lower temperature of 
spectrum A.2 is partially offset by a larger radius (Table[3}- It is 
remarkable that R in the models with reduced hydrogen and he- 
lium contents drops monotonically with time; in the basic model 
R increases between spectrum A.2 and spectrum A.3. 

To investigate whether it is possible to ascribe the change be- 
tween spectra A.2 and A.3 with a change in the blackbody and 
collisionally ionized emitters only, we fitted spectrum A.2 with 
an alternative model, in which all parameters related to the ab- 
sorption components are fixed to the values for SpectrumA.3, 
and only the parameters for the blackbody and for the collision- 
ally ionized plasma are allowed to change. The resulting fit is 
significantly worse (x 2 5000, as compared to the best basic fit 
X 1 — 4000). The^- 2 for this fit is dominated by the bad fit to the 
emission between 25-29 A, and around 16 A. 

To investigate whether the lower luminosity may be ex- 
plained by increased absorption due to increased columns of the 
ionized shells only, we fit spectrum A.2 with a model in which 
the parameters for the central emitters are fixed to the best values 
of the basic model for spectrum A.3, and only the columns and 
ionization parameters of the shells are allowed to vary. The best 
fit with this variant model has y 1 — 5500, significantly worse 
than the basic fit for spectrum A.2. 

To investigate whether the lower luminosity may be ex- 
plained by increased absorption due to an increased column of 
the neutral gas and/or an increased column of the dust, we fit 
spectrum A.2 with a model in which the columns of the neutral 
gas and of the dust are fitted, but all other parameters are fixed to 
the best-fit parameters of the basic model for spectrum A.3. This 
gives a bad fit, with x 1 — 8000. A model in which the central 
emitter is only partially absorbed by an increased neutral hydro- 
gen column does not fit the spectrum at all, due to the fact that 
the H I absorption is strongest at A > 24 A, a part of the spectrum 
that remains relatively bright. Thus, we can exclude absorption 
by cold clumpy material. 

Finally, we have checked whether an addition of a radiative 
recombination continuum (RRC) to the model of spectrum A.2 
improves the fit. The presence of this emission is suggested by 
emission features near 16.8 A and 18.4 A. We take into account 
the recombination of the most relevant ions of carbon, nitrogen 
and oxygen. The addition of RRC to the model does not lead to 
a significant improvement of the fit. 

4.3. Derived mass loss rates and model consistency 
Our fits to the spectra imply a mass loss rate: 



M w = \^p nr ^m H f c n H v = J47ifim H v- 



(2) 



where fi is the mean atomic weight in units of the hydrogen mass 
niH, f c is a factor allowing for dumpiness of the density distri- 
bution (f c < 1), and where Q. = 4n in the case of spherical 





layer 1 


layer 2 


layer 3 


M/(/ c £^)(M /yr) 

r/(/ cJ 8r BB ) 

n H X {f c pf 


1.3 x 10~ 5 

127 
2.2 x 10 18 


1.4 x 10- 4 

202 
1.2 x 10 19 


4.1 x 10- 2 
2.6xl0 7 
1.8 x 10" 



The parameters listed are valid for the standard model for spectrum A.3. 
If the layers are depleted in hydrogen and/or helium the actual mass loss 
rate would be lower. 

symmetry. The observed spectrum gives £ and v for each shell, 
almost independent of L. With the parameters of the standard 
model fit to spectrum A.3, we obtain the mass loss rates listed 
in Table [4] Note that these values are valid for one time interval 
only - i.e. the time of maximum X-ray flux, and thus may over- 
estimate the outflow-rate averaged over the outburst. From Eq.[2] 
we note in addition that the derived value of M w scales with the 
luminosity L of the central source, and thus will be smaller when 
the layers are hydrogen- and/or helium-poor (see Table 3), and 
also when the layers are spherically symmetric, in which case 
the xabs model overestimates the central luminosity. 

The depths of the absorption lines, combined with the fitted 
abundances, are converted to the hydrogen column density Nh = 
J nudr =s f c rinAr where AR is a typical length scale. Our model 
implies that this length scale is smaller than the distance to the 
central source, Ar = fir < r, and this may be used to obtain an 
upper limit to the distance r. 



Hi 

N H 



fcP 



fc n H Pr 

and from this a lower limit to the density 



<fc,P< 1) 



n H 



$r 2 L \f c fi 



(3) 



(4) 



Note that the definition and computation of £ involves the physi- 
cal density (i.e. not the space-averaged density), and therefore 
does not contain the dumpiness factor f. . The mean atomic 
weight fi is 1.42 for proto-Solar H-He abundances and 3.58 
for H-He abundances decreased by two orders of magnitude. 
Another interesting parameter is the time scale at which the 
gas in the outflow adjusts its ionization balance. This may be 
calculated from the ionization balance and ^-values of the fits. 
Through the SPEX tool recJime we determine the O vm recom- 
bination times for layers 1 and 2. They spread between a few 
(H-He rich slabs) to up-to-ten (H-He poor slabs) seconds, which 
means that the gas is responding instantaneously to the varia- 
tions of the photo-ionizing continuum of the source and that the 
observed features are linked to the variability of the source. 



5. Discussion 

From our fits it is clear that the central blackbody emission 
varies, during our first observation, rapidly, on a time scale of 
hours. The variations are likely exaggerated by the slab model: 
in full spherical symmetry, and also in the absence of hydro- 
gen and helium in the ionized absorption shells, the changes are 
still significant, but much less. Forty days into the nova out- 
burst, we may well be looking at shells whose abundances are 
affected by nuclear reactions. Indeed, the abundances that we 
derive for carbon, nitrogen and oxygen indicate that such is the 
case. This implies that the hydrogen and helium abundances of 
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absorption shells are much lower than solar, in comparison to the 
elements just mentioned. Even in the absence of spherical sym- 
metry, therefore, the variation of the central source is less pro- 
nounced than implied by our basic model, and probably more in 
agreement with the numbers found for the models in which the 
absorption shells contain no hydrogen and/or helium. 

The response time required for the ionization equilibrium of 
shells to respond to changes in the ionizing continuum is a few 
seconds, which means that the photo-ionization balance of the 
shells responds almost instantaneously to variations in the source 
flux and that its luminosity might have indeed changed. The 
mass loss rate in our standard model for spectrum A. 3 for the 
third shell in particular is rather high compared to the common 
limits 10~ 4 - 10~ 3 M o /yr found for other novae (|Kovetz et all 
119871: iKato & Hach isu 1989: ISmith & Owockil200fih . The de- 
rived mass loss rates for all three layers are lower in the models 
with reduced hydrogen and/or helium content of the absorbing 
shells. 

As discussed in Section 4.1, the blue-shifted absorption lines 
are caused because photons from the central source moving ini- 
tially in our direction are scattered into another direction. If the 
dimension of the shell is large with respect to the size of the 
white dwarf, and if the optical depth of the shell in the lines is 
not too large, this will lead to a P Cygni profile, as we see pho- 
tons initially moving in other directions and then scattered into 
our direction in emission, both blue- and red-shifted. The emis- 
sion and absorption should have equal strength. In the spectrum 
of V2491 Cyg this is clearly not the case. One explanation is 
that the shell is not large with respect to the white dwarf, so that 
a large fraction of the receding half is occulted. Similarly, a high 
optical depth in the shell may cause photons scattered from the 
receding half of the shell to be reabsorbed before escaping in our 
direction. And finally, a large asymmetry in the shell also may 
reduce the strength of the red-shifted emission. 

Carbon and sulfur are about an order of magnitude less abun- 
dant than nitrogen and oxygen, and elements as silicon, argon 
and calcium are much less abundant. This suggests that the white 
dwarf in V249 1 Cyg is a O-Ne rather than a C-O white dwarf. 
The strength of the Ne ix-x lines as compared to the C v-vi lines 
supports this conclusion, a s do observations in other wavelength 
regions dLvnch et al.ll2008t iHelton et al.l[2008t iNaik et al.ll2009l 
iMunarietal .11201 lh . 

The strongest feature in all X-ray spectra is the rest frame 
Oi K-edge around 23.0 A, which is well reproduced in the fits 
together with the corresponding Oi ls-2p absorption line at 
23.5 A. Other features of the cold CSM/ISM are the Fei L- 
edge around 17.5 A and the Ni l-2p absorption line at 3 1.2 A, 
both well fitted by our models. Neutral oxygen and nitrogen 
are highly overabundant, whereas iron is slightly overabundant, 
with respect to solar abundances. This suggests that part of the 
cold absorption originates in enriched circumstellar matter close 
to V2491Cyg. The reduction of the column and the decrease 
in temperature of the cold absorber between spectra A. 1 and B 
suggests that the CSM column is mostly from the current nova 
outburst, and still expanding, albeit slowly. Page et al. (2010a) 
suggest that the reduction of the column stops at 2.2xl0 25 rrT 2 
some time after our second observation. This final value then 
corresponds to the interstellar component. If the total O i column 
is equally divided between CSM and ISM, the ISM may have a 
near-solar abundance for oxygen, and sub-solar abundance for 
iron, implying that iron is depleted by dust. 

The optical flux defines the nova outburst, and thus clearly 
must come from the expanding nova envelope; the ultraviolet 



flux changes in tandem with the optical flux, and therefore also 
originates from the nova envelope. As can be seen from Fig. [5] 
our model fails to predict the correct level of the ultraviolet and 
optical fluxes, by several orders of magnitude. This indicates that 
our model applies only, if at all, to the X-ray part of the spec- 
trum, which must originate from a very different layer in the 
expanding envelope, presumably much closer to the white dwarf 
than the layers responsible for the optical/ultraviolet spectrum. 
This conclusion is confirmed by the observation that the marked 
changes in the X-ray flux during the first XMM-Newton obser- 
vation are not accompanied by changes in the optical and/or ul- 
traviolet fluxes. 

Because the absorption lines from the shells are shifted from 
the rest wavelengths by the high outflow velocities, and because 
hot atmospheres do not show deep absorption lines even at the 
rest wavelengths, the use of a blackbody spectrum as the central 
source is not as bad as might appear at first sight. As a first test 
we have tried a model atmosphere spectrum that replaces the 
blackbody as the central source. Following Ness et al. (2011) 
we ha ve used the best fitting TMAP model of iRauch & Werner! 
d2010h (007, with logg = 9 and T = 10 6 K). This provides an 
even worse fit with^ 2 = 15 with three shells, 18 with two shells, 
32 with just one shell and 59 without any photo-ionized absorber 
(instead of 7 from the BB fit) due to the fact that O vm / N vn line 
ratio (the most important in the spectrum) cannot be reproduced 
by the atmosphere model. It may be necessary to consider mod- 
els with different abundances and outflows which, however, are 
currently not available and we therefore refrain from a system- 
atic parameter study of the atmosphere model parameters. 

6. Conclusions and prospects 

Our spectral analysis of nova V2491 Cyg suggests that the ab- 
sorption by highly-ionized ions is due to non-monotonous, pos- 
sibly discrete, ejecta shells with different outflow velocities and 
ionization levels. Variations on time scales of hours occur both 
in the luminosity of the central source, and in the ionization level 
and columns of the absorption shells. We find that in our scenario 
of photo-ionization equilibrium the expanding gas is responding 
almost instantaneously to the variations in the source ionizing 
continuum. Our upper limits for the mass loss rate of each shell 
are in agreement with those estimated for other novae, especially 
if the shells are depleted in hydrogen. 

The values derived from the fits for the luminosity and ra- 
dius of the central source depend on the nature of the models 
that are applied, and thus must be considered uncertain. In the 
spherically symmetric case, however, the radius of the central 
emitter of X-rays decreases monotonically with time, and thus 
must be in the increasingly transparent expanding nova enve- 
lope. The values derived for the abundances in the ionized ab- 
sorption shells do not depend on the details of the models, and 
thus may be considered more secure. These abundances indicate 
that the white dwarf in V2491 Cyg is an O-Ne white dwarf. 

Improvements to our pilot models can be made in several 
ways. First, it would be useful to replace the central blackbody 
emitter with an appropriate white-dwarf atmosphere model. We 
have tested a model in which the TMAP model substitutes the 
blackbody, but still absorbed by the three shells. However, this 
model provides results even worse than those obtained with our 
standard fits with a blackbody continuum. In the future atmo- 
sphere models must take into account different abundances and 
outflows. 

Second, it will be necessary to iterate between the computa- 
tion of the ionization structure of the ionizing shells for a given 
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spectral energy distribution, and the computation of the spec- 
tral energy distribution and elemental abundances from the fit, 
for each of the four observed spectra separately. For this im- 
provement it is necessary also to determine the distance more 
reliably: a close distance implies a lower luminosity and thereby 
a lower ionization parameter of the shells. Finally, it would be 
useful to obtain constraints on the hydrogen contents of the hot 
ionized absorption shells, perhaps from ultraviolet observations. 
Ultraviolet observations may also help in determining the loca- 
tion of origin of the optical/ultraviolet flux. 
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